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Abstract 

We consider electron transport through a mobile island (i.e., a nanomechanical oscillator) which 
can accommodate one or two excess electrons and show that, in contrast to immobile islands, the 
Coulomb blockade peaks, associated with the first and second electrons entering the island, have 
different functional dependence on the nano-oscillator parameters when the island coupling to its 
leads is asymmetric. In particular, the conductance for the second electron (i.e., when the island 
is already charged) is greatly enhanced in comparison to the conductance of the first electron in 
the presence of an external electric field. We also analyze the temperature dependence of the two 
conduction peaks and show that these exhibit different functional behaviors. 

PACS numbers: 85.85. +j, 73.23.Hk 
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I. INTRODUCTION 



nano- 
1. In 



Electron transport in nanoelectromechanical systems (NEMS), such as suspendec 
beams, cantilevers, and nano-oscillators, is now attracting considerable attention 
shuttles, electrons can be carried by a single nano-particle or single molecule, which oscillates 
between two leads. This mechanical motion strongly modifies the lead-shuttle tunneling 
matrix elements, affecting the charge transfer. Theoretical j^, [], J] and experimental [5, Q 
studies of nanooscillators clearly demonstrated the influence of mechanical motion on their 
electrical properties. 

Previously, electron transport through a moving island was examined in the strong 
Coulomb-blockade regime, when the conducting level of the nano-oscillator can only be sin- 
gle populated, with higher-energy states being energetically inaccessible. Here, we demon- 
strate that a charged island behaves differently from an uncharged one; correspondingly, the 
possible double occupation of the conducting level leads to a situation where the Coulomb 
blockade peaks, associated with the first and second electrons transferred through the island, 
have different dependencies on the nano-oscillator parameters. Moreover, we show that the 
double occupation leads to a conductance enhancement for the second electron entering the 
island. To achieve that, we apply a previously-developed approach [3, Q] which makes it 
possible to examine the case of finite on-site Coulomb interaction. 

It should be noted that the moving island studied in this work can be considered as a 
shuttle because of its actual function: shuttling. When the island moves closer to the left 
lead, increasing that matrix element, it loads an electron. Then, the island moves closer 
to the right lead and unloads the electron. Therefore, this describes the operation of an 
electron shuttle. We have opted to use the term "mobile island" for this "electron shuttle" 
because in the nanomechanical community oscillators exhibiting an instability are called 
"shuttles". Still, functionally, the moving island described here is effectively a shuttle. 

Usually, a nano-oscillator is considered to be placed symmetrically between the leads. 
However, recently several works discussed the situation when there is an asymmetry in 
;he lead-oscillator coupling produced either by the difference in the tunnel matrix elements 
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10j or by the spatial shift of the equilibrium oscillator position In the latter case, 



it was theoretically proposed that, if the island is closer to one lead than to the other, the 
current through the structure depends exponentially on this spatial shift (with the tunnelling 
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length A), because the overlap integral of the electron wave functions in the island and in 
the leads, involved in the tunnel matrix elements, exponentially decreases with distance. 
In the model of Ref. 11|, this small displacement, shifting the island close to one of the 



leads, was produced by the large magnetic field gradient acting on the spin of the nitrogen 
or phosphorus impurity incorporated into their model of a Ceo shuttle. Here we show that 
such displacement can be achieved naturally in the island without impurities with an excess 
electron in an electric field produced by the source-drain voltage or an external capacitor (see 
Fig. 1). Moreover, this kind of spatial asymmetry can be associated with the Jahn- Teller 
effect: when an orbital state of an ion is degenerate for symmetry reasons, the ligands will 
experience forces driving the system to a lower-symmetry configuration, lowering its energy. 
Consequently, the ligand position between the two ions is not symmetric and changes with 



the electron transfer 
atoms in manganites 



rom one ion to the other. Oscillations of suc h lig ands, either as oxygen 



or rare-earth atoms in filled skutterudites 



13], were analyzed jointly 



with the Jahn- Teller effect. However, the tunnelling length was assumed to be infinite in 
Refs. Q,0| and the dependence of the tunnel matrix elements on the oscillator position 
was not taken into account. Here we consider these effects and find a remarkably rich 
behavior of the conductance of nano-oscillators, if the matrix elements have an asymmetry 
as in Refs. 0,Q. 

The present paper is organized as follows. Section II introduces the pertinent Hamiltonian 
including all interactions. The equations of motion for the electron creation/annihilation 
operators are derived in Sec. III. The equations for the electron populations and the popu- 
lations correlator are derived and solved in Sec. IV. In Sec. V we obtain explicit expressions 
for the lead-to- lead current and discuss the dependence of the conductance on the system 
parameters. The conclusions of this work are presented in Sec. VII. 



II. FORMULATION 



To examine electron transport through a moving island, we assume that the island has 
a single spatial state which can be populated by two electrons having opposing spin projec- 
tions, a and a, with finite on-site Coulomb interaction (Uq ^ oo). It should be noted that 
here we consider the situation where the coupling of the nano-oscillator to the leads is weak, 
so the Kondo-like correlations are not important. The Hamiltonian of this system is given 
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by (a — L,R for left, right; a = 1, 2 for spin up, down; a = 2, 1 ) 

H = ^ E a a t a a + U o a t a a a t a a + Ekaa °k<*° Ckaa + Hosc + ^ tun ' ^ 

a kaa 

where a+ (a a ) are the creation (annihilation) operators for the electrons in the island and 
c taa ( c kaa) are the creation (annihilation) operators with wavevector k in the a-lead. The 
tunneling term, 

#tun = - ^ Tka Wa ^> °kaa a a +H.C, (2) 
kaa 

has tunneling amplitudes depending explicitly on the position x of the island as 

w a (x) = exp 

with the tunneling lengths A^ = —A and Ar = A for the left and right leads, respectively. 
The Hamiltonian of the nanomechanical oscillator also contains the interaction between the 
charge stored in the oscillator and an effective electric field £, as 

i*» = £ + ^ (3) 

a 

This field S can be produced by the voltage applied to the leads, by the Jahn- Teller effect, 
or even by an independently controlled electric field, if the structure is placed inside an 
external capacitor. Here, N a = a+ a a is the electron population operator, and M and ojq are 
the effective mass and the resonant frequency of the nano-oscillator, respectively. 

After the unitary transformation U = exp{— ip ^2 a X£ N a }, where X£ = eS /(Muq), we 
obtain the usual expression for the oscillator Hamiltonian, H osc = p 2 /2M + MuJqX 2 /2, and 
the modified electron operators, 

a' a = U + a a U = e- ipx °a a , (4) 

and tunnel matrix elements, 

w' a (x) = w a (x + ^x e N a ,). (5) 

a' 

Using the properties of Fermi operators (N 2 = N a , a a N a = a a , N a a a = 0, ...), we obtain 

w a (x + x £ (N a + N & )) = w aa (x) + [w a (x + x £ ) - w a (x)]N a 
+[w a (x + x e ) - w a (x)]N„ + [w a (x + 2x £ ) - 2w a (x + x e ) + w a (x)]N a N & , (6) 
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and 

a a w a (x + x £ (N a + N 9 )) = a a w a (x + x £ ) + A a [w a (x + 2x £ ) - w a (x + x £ )\. (7) 

Here, we have introduced the Fermi operator A a = N„a a . Accordingly, the tunneling term 
has a form 

#tun = - T ka c\ aa B aa - H.c, (8) 

kaa 

where B aa is the Fermi operator given by 

B aa = a a (u a + v a Na), (9) 

with the first term responsible for the electron tunneling from the unoccupied island and 
the second one describing the electron transfer through the double-populated level. Here, 

u a = u™ n p mn and v a = v™ n p mn , (10) 

where 

Pmn = \m)(n\ (m,n = 0,1,..) (11) 

are the eigenstates of the mechanical oscillator Hamiltonian, and the matrix elements of the 
tunneling amplitudes are given by 

w™ n = (m| exp f — J ex-p(—ipx £ ) \n) (12) 



and 

«r = < n [exp ( f ) ~ l] ■ (13) 

n 

Equations (12,13) can be considered as a generalization of the Frank-Condon factors [141 ] 
accounting for the overlap integral of the vibrational states before and after the transition. 
It is evident that the Frank-Condon factors are different for the first and second electrons 
entering the island because the center of the oscillations is shifted in the case of the charged 
island (see Fig. 1). It should be emphasized that by introducing the operators B aa we are 
able to derive the equations of motion analytically, without the use of the Hartree-Fock 
approximation, assuming only a weak lead-island tunnelling coupling. From a general point 
of view, the method presented here is equivalent to the master equation approach. 
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III. EQUATIONS OF MOTION 

Equations of motion for the island electron operators obtained from the Hamiltonian, 
Eq. (1) are given by 



ia a = E a a a + U A a - ^T fc * QCT (u+ CT + v+ a N„)c kaa + 

kaa 

T ka*4 a ^a a aaV a v + Y T kaa V aa a t a v C kaa (14) 



kaa kaa 



and 



iA a = (E a + U )a a - J2 T kaa«a + V aa)^^ C kaa + 

kctcr 

5^ T k aa4 a ^a a u a5 + ^ T^(ul„ + v^)ata a c ka& . (15) 

kaa kaa 

Accordingly, equations for the ensemble averaged island populations can be written as 

^ = ~i T ^MaaBaa) + H. C. (16) 
kaa 

and 

d =^Y1 T kaa(4 aa A aa ) ( (U aa ) + (V aa )) - i ^ T ka9 ( C £ aff A a& ) ( (u a& ) + (v a& )) + H. C. 

kaa kaa 

(17) 

The equation of motion for the electron operators in the leads are given by 

i&kaa E kacr C ka<J T kacr B acr . (18) 

In the case of weak lead-island tunnel coupling, the solution of this equation can be repre- 
sented as 

Ckaa(t)=C k °l(t)-T kacT J dhgl^it^Baaih), (19) 

where c k °^ a (t) is the unperturbed electron operator and g ka(7 (t,ti) * s the retarded Green 
function of the electrons in the leads, given by 

9l a aitM) = = -ie- iE ^9(r), (20) 

where [...,...] + is the anticommutator, r = t — t 1 , and 6{r) is the unit step function. It 
should be emphasized that the non-Markovian dynamics involved in Eq. (19) allows us to 
reveal manifestations of the oscillatory mechanical motion during the tunneling events. 
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IV. ELECTRON POPULATIONS AND POPULATIONS CORRELATOR 



A. Free-evolution approximation 

To determine electron populations in the island and the correlator of the populations 
having different spin projections, we substitute Eq. (19) into Eqs. (16,17). The correlators 
of the type (cf^ (t)B aa (t)) can be rewritten using the formula 

E T ^<^(*)^(o> = -EiW f ^i(41 + w41(*i)>([^w,^(^)] + >. (21) 

kaa kaa 

To decouple the correlators for the electron operators in the island, we use the approximation 
of their free evolution, which is valid in the case of weak lead-island tunneling. The free 
evolutions of the operators a a , A a , and B aa are given by 

a a (t) = e-^Mfx) - (1 - e-™° r )Mti)], (22) 
A a (t) = e-^+^A^U), (23) 

and 

B aa (t) = e-^Mh) - Mti)Mt) + e-^+^A^iXxw + v aa )(t). (24) 
Accordingly, 

(B aa {t)BZ a {h)) = e-^K CT (t) M + (T (^i))(l - N a - N- a + N a N- a ) + 

e- i{EM)T ((u aa (t) +v aa (t))(ut a (t 1 ) +v+ a (h)))(N & -N a N 9 ), (25) 

and 

(B+ (*i)£~(f)> = e- iE ° T {u+ a (ti)u a a(t))(N a - N a N & ) + 

e~ i{EM)T ({utiti) +vL(ti))(u*a(t) +v aa {t)))(N N 9 ). (26) 

The free-evolution approximation can also be used to calculate the correlators of the me- 
chanical operators. Using p mn (t) = e 1 "" 171 ^ 1 " 1 p mn (ti) , we obtain: 

(««r(*)«^(*l)>=El U -| 2e * rtm " T <An>. 
mn 

(^waoi = E \<:\ 2e ~ iulmnT (pn), (27) 

mn 

where (p n ) = (p nn ) is the steady-state distribution of the mechanical degrees of freedom and 
w ra « = e m - e n = uj (m - n). 



B. Electron occupations 



In the absence of an external magnetic field, the averaged electron populations, (Aq) and 
(A 2 ), should be equal: (Ni) = (N 2 ) = (A). As a result, we obtain the following equations 
for the averaged electron occupation (A) and for the correlation function of the populations 
with opposite spin projections, (NiN 2 ), 

Vl (N 1 N 2 )=r ]2 (N), V3(N)+ V 4(N 1 N 2 )= Vo , (28) 

having the simple solutions 

(AQ = 770771 , ( Nl N 2 ) = W2 . (29) 

V1V3 + V2V4 vm + V2V4 

We introduce the following coefficients 

a mn 

Vl = J2J2 r a K n \ 2 e 2X£/Xa [(Pn) + (Pm ~ Pn)f a (E + U - U mn )], 
a mn 

V2 = J2Y, T «\< n \ 2e2X£/Xa ^f^ E * + U °~ U rnn), 
a mn 

V3 = Yl Yl r al U a"| 2 [(Pn> + (Pm ~ Pn) f a (E - U mn ) + 
a mn 

(p m )fa(E - U mn ) - 6 2X£/Xa (p m )f a (E + U ~ UJ mn )}, 

?74 = YY Tc <\ U ™ n \ 2 ^P n -pm)fa(E ~ UJmn) ~ 
a mn 

e 2x£/x «( Pn - Pm)f a (E + U - uj mn ) + (e 2x£ / x « - l)(p n }}. (30) 

Here, f a (E) are the electron Fermi distribution functions in the corresponding lead and (...) 
means ensemble averaging. In the wide-band limit, we can introduce the tunnel rate as 

T a(T = 27r ^ \T kaa \ 2 5(u - E kaa ). (31) 
k 

In this work, we examine the case of a very small source-drain voltage applied to the system, 
so the density matrix of the mechanical oscillator has the equilibrium form (ks = 1) 

Pm = e -^om/T ^ _ e -^o/T) _ (32) 

We plot the solutions, Eq. (29), as well as the second cumulant, 

K N = (A^) - (iV> 2 , (33) 



in Fig. 2 as functions of the separation between the energy of the island level E a and the 
equilibrium chemical potential fi of the leads. The following set of parameters, associated 
with C 60 , was chosen j^j]: the charging energy, U = 270 meV, the fundamental frequency, 
TujJq = 5 meV, and the fundamental uncertainty of the oscillator position 



The magnetic field is taken to be zero (so E a = Eg), T = 77 K, eV = 0.5 meV, and 
A = 4 pm. It is evident from Fig. 2 that when the electron energy level on the island 
becomes smaller than fi (modulo thermal broadening), the island is single-populated and, 
when the energy separation between E a and fi is larger than the charging energy, the island 
is double-populated, as expected. It should be emphasized that although the ensemble 
averaged values of both electron populations are nonzero in the case of single occupation, 
the population correlator is zero, meaning that the electron having only one of the spin 
projections can be found in the specific sample. This Pauli repulsion also manifests itself in 
the negative value of the cumulant K N in the single-occupation regime. It should be also 
noted that the functional dependencies of Fig. 2 do not depend on the value of xg and the 
asymmetry of the couplings to the left and right leads. 

V. ELECTRON CURRENT AND CONDUCTANCE 

The current flow of electrons having a-projection of the spin from the a-lead can be 
defined as 



r = ^Jh/2Muj = 3.8 pm. 




(34) 



k 



k 



Using the same approximations as in the previous section, we obtain 



I a eV a ^ 



u. 



mn 1 2 



a 



{(1 - f a (E - u mn ))(N - NiNzXpn) 



rnn 



+e 



-f a (E - u mn )(l - 2N + iVxiV 2 )(p m ) 
C£ / A "[(l - f a (E + U - w mn ))(JViJV 2 >(p n > 
-UEo + Uo-u^iN-Nt^ipm)}} . 



(35) 



The associated conductance, 



(36) 
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is presented in Fig. 3 as a function of (Eq — /i) and xs, using the same parameters as in 
Fig. 2 with coupling constants KFl = 0.1 meV and hT^ = 0.002 meV and temperatures 
(a) T = 4 K and (b) T = 77 K. The projections of the three-dimensional plots unto both 
the U G versus (E a — fi)" and "G versus x?" planes are shown in Fig. 3(c). One can see 
from Fig. 3 that the magnitudes of the conductance peaks, associated with the first and 
second electrons entering the island, are only equal to each other for xs = (conventional 
Coulomb blockade case). Moreover, the conductance through the charged island is drastically 
enhanced at positive moderate values of x&. It should be noted that the electric field- induced 
shift would not produce a conductance enhancement for the immobile island because the 
exponential increase of the tunnel matrix element between the island and one of the leads 
is compensated by the same exponential decrease of the tunnel matrix element coupling to 
the other lead. However, for the mobile island, these matrix elements are averaged over the 
island oscillatory motion and the shift is not cancelled out. This is even more pronounced 
for the charged island where the center of the oscillations is already shifted by the presence 
of the first electron. Formally, the account of the oscillatory mechanical motion during the 
tunneling events becomes possible due to the non-Markovian character of the equation of 
motion. The dependence of the conductance peaks on x? has a Gaussian form (coming 
from the Frank-Condon factors) with the centers shifted to two different positive values 
of X£. With increasing temperature, the conductance peaks become broader and the shift 
is increased, as seen in Fig. 3. This shift can be attributed to the phono n-blockade effect 
discussed in Refs. It should be noted that the bias-voltage independent displacement 

xs can be created, for example, by nearby charge impurities, image charges, device geometry, 
etc. Q. 

We also examine the temperature dependence of the conductance peak magnitude. For 
the immobile island, one can expect either no temperature dependence, in the case of 
quantum-mechanical tunneling, or thermal-activation dependence, in the case of over-the- 
barrier hopping transport. However, deviations from such behaviors were observed both in 
transport through single molecules 15] and in the resistivity of manganites [la] • Theoreti- 
cally, it was shown that either the mechanical motion of the nanoconductor or coupling to 
the quantized thermal modes can lead to exotic types of temperature dependence. 

These can be seen in Fig. 4 for various values of A, at x$ — 0. It should be noted that the 
curves are identical for both peaks in this case. For nonzero X£, the temperature dependence 
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becomes even more complicated for small A, as can be seen in Fig. 5(a) and 5(b) for the 
charged and uncharged shuttle, respectively, because of the temperature-induced shift of 
the peak position (see Fig. 3). It should be emphasized that the functional dependencies 
at small A are very different for the two conduction peaks, with the peak for the second 
electron being almost an order of magnitude larger. 

The current l aa through the shuttle is extremely sensitive to the value of A, because it 
appears in several exponents of Eq. (35). It is evident from Fig. 4 that the smaller A is, the 
larger the conductance of the system becomes. Therefore, the quality of the leads plays a 
more important role in the electrical properties of nano-oscillators than in most standard 
electronic devices. 

VI. CONCLUSIONS 

In conclusion, we have examined electron transport through a mobile island which can 
contain one or two electrons. We have derived the equations of motion for the electron 
creation/annihilation operators and have been able to evaluate them (without using the 
Hartree-Fock approximation) by introducing a complex Fermi operator B ar7 , Eq. (9). Based 
on this microscopic approach, the equations for the island populations and the correlator of 
the populations have been derived and solved. They are involved in the expression for the 
electron current through the structure, also obtained microscopically. We have shown that 
in the presence of an external electric field (produced either by the voltage applied to the 
system, by the Jahn- Teller effect in the molecular junctions, or by an external capacitor), 
and an asymmetry in the coupling of the island to the leads, the conductance of the second 
electron entering the charged island is greatly enhanced. The temperature dependence of 
the conductance has been also discussed. 

This work was supported in part by the National Security Agency, Laboratory of Physical 
Sciences, Army Research Office, JSPS CTC Program, and National Science Foundation grant 
No. EIA-0130383. 
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(a) (b) 




| ►! X£ ,< 

FIG. 1: (Color online) Schematic diagram of electron transport through (a) uncharged and (b) 
charged nano-oscillator (electron charge e is negative). 
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FIG. 2: (Color online) Electron populations (for spin up, a = 1, in red, for spin down, a = 2, 
in dashed black), populations correlator (dotted-dashed intermediate purple curve), and cumulant 
(green dotted curve at the bottom) as functions of the energy of the nano-oscillator electron state. 
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FIG. 3: (Color online) Conductance of the nano-oscillator as a function of the state energy (E a — [/,) 
and the oscillator shift xg in an external electric field £, with a tunneling length A = 4 pm, hTz, = 
0.1 meV and hF R = 0.002 meV for (a) T = 4 K and (b) T = 77 K. (c) Left: projections to the "G 
versus (E a —fj,)" plane. Right: projections to the U G versus xg v plane. The dashed red (continuous 
blue) peak pi (P2) denotes the conductance peak at E a = fi (E a = \i — Uq). When xg = 0, pi and 
P2 have an equal conductance, corresponding to the usual Coulumb blockade results. 
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G ~ constant 



100 200 300 

T (K) 



FIG. 4: (Color online) Temperature dependence of the conductance G(T) for various tunneling 
lengths A and for an oscillator shift = 0. The yellow (shadowed) ellipse is the low-temperature 



regime where G ~ T 
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FIG. 5: (Color online) Temperature dependence of the conductance G(T) (for various tunneling 
lengths A and for an oscillator shift xg = 5 pm) for (a) pi at E a = fi and (b) p2 at E a = \x — Uq. 
The conductance can vary orders of magnitude for small changes in A. 
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